Regularisation as a quantised low-pass filter 
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A divergence-free approach to relativistic quantum electrodynamics based on regularisation of 
equations of quantum mechanics is discussed. This approach is shown to be exactly equivalent to 
the conventional Feynman-Dyson renormalisation techniques. 
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Introduction. — It is appropriate to start this letter with 
a word of caution. This is not a solution to the problem of 
infinities in quantum field theory. My goal is to construct 
a computational scheme with the standard perturbation 
approach leading directly to regularised diagram series, so 
that infinities never occur (for all standard definitions see 
Refs. 0,0)- This scheme may be regarded as a construc- 
tive implementation of Bogoliubov's point that equations 
of quantum field theory are mathematically "underde- 
fined" and hence contain some freedom. This freedom 
is implemented by defining "observable" fields as mix- 
tures of the physical and fictitious ones. The interaction 
is then introduced in such a way that divergences never 
occur. This yields a perturbation expansion which coin- 
cides exactly with the Pauli-Villars-regularised diagram 
series in the conventional approach [l[ , with the fictitious 
masses becoming regularisation masses. The standard 
renormalisation procedure then "traps" the freedom in- 
troduced into the equations of motion in observable mass 
and charge of the electron. Unlike in the conventional 
Feynman-Dyson approach, all expressions remain math- 
ematically meaningful at all stages of calculations. 

Response and regularisation. — In detail, this connec- 
tion will be discussed elsewhere [3JI, see also 01 . The key 
idea may in fact be demonstrated for the oscillator with 
a variable frequency, with the Hamiltonian, 

H v (t) = [fi -«(<)] ala . (1) 

Here, ao,aJ is the standard creation/annihilation pair in 
the Scrodinger picture (SP), [ao,&o] = 1, Oo is the os- 
cillator frequency, and v(t) is a c-number function. We 
use units where h = c = 1. Operators in the interac- 
tion picture (IP) will be signalled by the presence of the 
time argument: aj(t), ao(t), and those in the Heisenberg 
picture (HP) by uppercase letters: A^t), A (t). These 
conventions also apply to the field operators below. 

Feynman's diagram technique [l| for this toy system 
comprises two graphical elements, 
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where (cf. Fig. [T]) 
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FIG. 1: Integration contours in the complex plane of energy. 
The closed contours encircle all poles of integrands (shown 
schematically as crosses). Poles on the negative real axis ap- 
pear in relativistic models. For details see e.g. Schweber's 
textbook pj|. 



The only divergent diagram Q in this technique is the 
contribution to the vacuum phase, with v = J dtv(t), 
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The logarithmic divergence here reflects the fact that 
Ap'°^(0) is not defined. One regularisation suffices to 
cancel the divergence: the regularised propagator, 
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is continuous at t = 0, and a1 0,1 '(0) 



0. 



The regularised propagator and loop can be expressed 
as legitimate diagrams in another diagram technique, 
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This technique comprises four graphical elements, 
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= O x , 
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(t-f), 
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v(t), (8) 



where Ap 1 ' 1 ^ is given by Eq. ([3]) with Oo — > Oi. Numbers 
on lines and vertices must match. Such diagram tech- 
nique follows from the time-dependent non-Hermitian 
Hamiltonian, (see endnote Q) 

H v {t) = f7 &( 0) a (0) + OiS^a^) 

-OiS^fyo) -v(t)S| 0) a(i), (9) 
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where a(o),&(i) and SL^St^ are dual operator sets, 



'LK; 



L,K = 0,1. 



(10) 
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L a (£)' V)- 
Heisenberg equations related to ([9]), 

(no-ift)i{o)(t) =«(i)i ( i)(t), 

(Q 1 -i5t)i(i)(t)=nii (0 )(t), 

implement a natural physical concept: a signal circulat- 
ing in a feedback loop, with regularisation related to the 
infinitesimal delay in the loop. 

Diagrams (J7|) , ((8]) imply that the first two terms in © 
are regarded as free Hamiltonian and the last two as in- 
teraction. What if we consider only the last term as inter- 
action and the sum of the first three as free Hamiltonian? 
To understand what such free Hamiltonian means, we 
replace the interaction term by the source term s(t)b\ Q y 

Equations (|11|) then hold with v(t)Am(t) — > s(t). Using 
them we find the linear response function of the system, 

(cf. Fig. EE! 
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It is instructive to compare this to Kubo's formula [B| 
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Quantum averaging present in Kubo's expression is omit- 
ted here because both sides are c-numbers anyway; the 
same applies to Eq. (TT2"|) . Manipulating complex con- 
tours in the conventional way [l| we recover the two-time 
commutator, (cf. Fig. [1]) 



[«(!)(*), &f )(*')] = 
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-ic(t-t') fix 
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Using the formula, (with a = 1/ ^/l — f2o/Qi) 
Hi cr 2 CT 2 



(Slo — o;)(f2i — u>) SIq 



(14) 
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we can express the dual operators in terms of the creation 
and annihilation operators of two harmonic oscillators, 



0(1) 



at 



= a(a a + ai), 
a(&l - a] 



ao, 



(16) 



i, nlt) 5} m = aj empha- 



The alternative notation S(x) = a ut, "( ) 
sises that these operators correspond to a specially chosen 
pair of input/output modes of the system. Unlike for a 
randomly chosen pair, Eq. (|12|) is regular: it lacks the 
step-function contribution. Using p6[) it is straightfor- 
ward to show that the free part in ([9]) is simply a sum of 
two oscillator Hamiltonians, 



H v (t)\ v=Q = iloalao + ttia\ai. 



(17) 



Toy model summary. — Regularisation is a result of (a) 
introducing fictitious degrees of freedom, (b) changing 
variables in the free Hamiltonian, and (c) modifying the 
interaction Hamiltonian by replacing a single physical 
field with a pair of fields with engineered commutational 
properties. In turn, these are inherently related to re- 
sponse properties of the system. Physically, the emer- 
gent regularised system is a quantised low-pass frequency 
filter with feedback. 

Regularisation algebra. — Let Ma, • • • , Mn (N > 0) be 
some real constants (all different). Generalising Eq. (TT5")) 

we define a set of complex functions G\ K ' L \z), (by de- 
fault all indices vary from to N) 



G\ K - L \z) 



M^l\i =K (^^/Mir\K<L, 



0, K>L. 
and real constants ek = ±1 and ok such that 



(18) 



N 
K=0 



Gf^(z) 



CO- 



<.M N , 



For the typical case Mq <C Mi <C 

e K = {-l) K , K = 0,- 
ct «°-i«1, < \a K \ <C 1, K = 2,- - ,7V. 

Of use to us will be the recursion relations (K < L) 



(20) 



[(M L ~z)/M L ]G^ L) ( 



[{M K -z)/M K+1 ]G\ K ' L \z) = G\ 



i(K+1,L) 



(z). 



(21) 



Dirac's filter. — Following the pattern of Eqs. ([TT]) . we 
wish to replace the equation for the electron operator in 
spinor quantum electrodynamics by the interacting Dirac 
low-pass filter, (with K = 1, ■ ■ ■ , N) 



(M - »0)*( O )(aO = e4(x)* {N) (x), 
(M fc - i0 {K) (x) = M K ${K-i){x). 



(22) 



Here, ^(^(x) are linear combinations of conventional 
spinor field operators *&i{x) with masses M;. Feyn- 
man's slash notation is used for a scalar product of a 
4- vector with Dirac's matrices |l[, 9 = = j^d/dx 11 , 
x = {t, r}, and 4 = 7 m ^a" wnere A\i is the 4- vector elec- 
tromagnetic potential. We proceed as for the toy model 
<j9j) , by firstly changing variables in the free Hamiltonian, 
and then modifying the interaction. The free Dirac low- 
pass filter is constructed from N+l free Dirac fields ipiix) 
with masses Mf. 

(Mi-i^i{x) =0, Z = (),•■• ,N. (23) 

Up to a factor 7 , these follow as Heisenberg equations, 

idMx) = $(z),#$(aO,^(a:);M,)]. (24) 
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where H(ip(x), i/)(x); M) is the Hamiltonian of free 
Dirac's field [l|], and ipi(x) = ijjj(x)j a is Dirac's conju- 
gate of ipi(x). The dual operator sets comprising the 
filter are given by recursion relations, (K = !,••• , N) 



N 



1=0 



N 

<f(o)(x) = $in{x) = ^2 tiai4>i(x), 



1=0 



M K <i>{K){x) = <P{k-i)(x)(M. 



K-l 



(25) 
(26) 
(27) 
(28) 



We now prove the all-important formula for two-time an- 
ticommutators, extending Eq. (|14p to fermions: 



c 



d A p 

H2-K) 4 



(29) 



where p 1 = 7 M p M , p = {po,p}- For K = N, L = 
Eq. (|29p follows from (JT9J) and from the complex-contour 
expression for the unregularised anticommutator [l[; for 
other K, L it is verified by making use of Eqs. (|2"Tj) . (12"6")) 
and (|28p . For t = t' it reduces to a counterpart of (TIT)]) , 



$ { K)(x),0 (L) (x')] + \ t = t , = 5 KL 5^{r-r')f. (30) 

The fields ip(K) (x) hence form a complete set, in the sense 
that the linear transformation between ipi (x) and ip/K) ( x ) 
is invertible. Using this we can extend equivalence of the 
Hamiltonians © and Q17p to Dirac's fermions: 



JV 



H(Mx),Mx); Mi) = H(0 ( o)(x)4lo)(n:);Mi) 

1=0 

N 

+ [^(kK)(x),4>(K)(x);Mi) 

K=l 

- M K0\ K) (x)i>(K-l) (x) (31) 

To start with, notice that each successive recursion in 
(121))) shortens the initial sum of operators ([23)1 by one 
addend, so that ip(o)(x) oc ^o^): and 



(M o -i0)$ {o) (x) = O. 



(32) 



Equations (|2"6")) . (|32p may be used as N + 1 equations of 
motion for the filter in place of (f2"3")l ; assuming otherwise 
contradicts the fact that ip(K)(t) form a complete set. 
Furthermore, the RHS of fl3T]) results Q in Eqs. ([H, 
(1521 . while the LHS— in Eqs. flU]). This is only possible 
if the RHS and the LHS of dSTJ) coincide. The "free Dirac 
low-pass filter" has thus been successfully constructed. 
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FIG. 2: Typical divergent sub-diagrams regularised by the 
"Dirac filter." 



Regularised quantum electrodynamics as quantised 
filter. — To close the "feedback loop," the Dirac filter 
must be supplemented by the modified electromagnetic 
interaction. The full Hamiltonian then reads (in IP) 

4eg(t) = H cl {t) + fl ph (t) + flint (t) (33) 
flmt(t) = -e / d 3 r^ in (x)4(x)tP out (x). (34) 



Here, H e \ is given by the RHS of (f3"Tj) . and H p h(t) is the 
standard Hamiltonian of the free electromagnetic field [H . 
Equations (|2"2"|) obviously follow from (f3"4")l . The equation 
for the electromagnetic field operator may be derived in 
a conventional way (with non-Hermitian current). 

Typical divergent diagrams regularised by the Dirac 
filter are shown in Fig. [2j where the graphical notation 
is an obvious generalisation of Eqs. ([7]), ©, cf. also end- 
notes @, [1] and Eqs. ([38]), (J39J) below. The electromag- 
netic propagator [l| remains unregularised. It is easy to 
see that N = 2 suffices to regularise all physically rel- 
evant diagrams, while N = 4 guaranties convergence of 
all diagrams, including the vacuum-phase ones and those 
subject to the Farry theorem Indeed, in any loop 
the number of electromagnetic propagators cannot ex- 
ceed the number of electronic ones. Hence the worst pos- 
sible divergence is exhibited by the "tadpole" (Fig. [2^); 
it is regularised with N = 4. The worst-case scenarios for 
physical diagrams are the self-mass of the electron (Fig. 
[2p) and the vacuum polarisation (Fig.[2j;); they become 
convergent with N = 2. The vertex correction (Fig. [2ji) 
is regular already with N = 1. 

Pseudo-adjoint and Keldysh techniques. — In Eqs. (|22p . 



fields 4 r (o)(a;), • • • , 4'(jv-i)( a; ) ma Y be eliminated result- 
ing in a closed regularised equation for a single field 
%) = $(N)(x): 



(M„-i0)(l-»WMi 



(1 - ip/M N )*{t) 

= e4(x)if(x) 



(35) 



This corresponds to eliminating intermediate vertices ([7]) 
in all propagators. The RHS in (TJ5*|) is the usual self- 
action source. Its action is infinitesimally delayed by the 
low-pass filter thus cancelling self-action of fermions at 
infinitesimally small times and distances. 

Intermediate fields can also be eliminated from pertur- 
bation calculations. It is convenient to introduce pseudo- 
conjugation of operators denoted by and the associated 
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indefinite scalar product denoted as (-||-), as 

X* = XX% ($ip 2 > = (*i|2"|$2>, (36) 

where |$i), ($2) and X are arbitrary state vectors 
and operator in Fock space, I = e^e" 1 is the 

sign operator, and are number operators for parti- 
cles with mass M;. All usual properties hold: (X+p = 
1, (aii + &1 2 )* = a*l| + b*xl (XiX 2 )i = X\X{, 
(<I>i||A > |<i>2)* = (<& 2 ||A >:,: |<i>i). A diagonal matrix element 
($||X*X|$)* = ($|i"tJX|$) is real but not bound to be 
positive. For brevity, we shall talk about Z-conjugation 
etc. The sign operator X does not commute with the 
interaction, so its use is confined to IP. 

By construction, yl)- m {x) is Dirac's Z-conjugate (de- 
noted by tilde) of i/) ut(x) = ^(x), 



N 



N 



1=0 



1=0 



The interaction Hamiltonian in IP is Z-Hermitian, and 
the S-matrix in the IP is X- unitary, (cf. also endnote Q) 



S = T + exp 
S = T- exp 



ie J d i x^{x) J^(x)'4>{x) 
— ie / d xtp(x)fi(x)iji(x 



(38) 



s- 1 . 



Wick's theorem "does not care" whether ip(x),tp(x) are 
genuine conjugates or X-conjugates, hence conventional 
perturbation calculations apply with the regularised con- 
traction, 



i(0|T+V>(x)^(x')|0) 



d 4 p 



leading directly to regularised diagram series, cf. Fig. [2] 
The two additional contractions characteristic of Perel- 
Keldysh's diagram approach Q, 

-i(0\i>(x)i>(x')\0), i(0\i>(x')<ip(x)\0), (40) 

are also given by regularised expressions, namely, by Eq. 
(1591) with C F -> C (±) (Fig. [J). A nontrivial problem is 
however posed by quantum averaging. For example, us- 
ing conventional perturbation techniques [1,0] we obtain, 

g(x,x';x",x'") 

= TrpT^ out (x)& m (x') T + $ oat (x")& m (x'") 

= TrpT_§ip(x)4(x')T + S4>(x")i>(x"'). (41) 

Here, p is the Heisenberg density matrix, or, which is 
the same, the initial state of the system. However, a 



consistent generalisation of the conventional perturbative 
expression implies also a repalacement of the quantum 
averaging by I-averaging, ((X)) — TrpIX, resulting in 

g{x,x';x",x'") 

= ((tM(x)^(x')T + S^(x'')^(x'''))). (42) 



This equations follows from ([41]) if we assume that the 
initial state contains only physical particles of mass Mq, 
so that p = pi. In turn, this assumption is consistent 
with the standard one P, Q , namely, that the interaction 
is switched on adiabatically, and that in the remote past 
the system is free. 

It should not be overlooked that we have just postu- 
lated an explicit distinction between the past and the 
future. This raises a question as to what extent regular- 
isations are consistent with reversibility. 

In conclusion, a divergence-free computational scheme 
has been constructed in the framework of relativistic 
quantum electrodynamics, and its equivalence to the 
conventional Feynman-Dyson renormalisation techniques 
has been demonstrated. 

The author is grateful to R. Glauber for triggering his 
interest in the topic, and to M. Efremov and Yu. Lozovik 
for enlightening discussions. 
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